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Abstract. The morphologies of detected jets in X-ray binaries are almost as diverse as their 
number. This is due to different jet properties and ambient media that these jets encounter. 
It is important to understand the physics of these objects and to obtain information about 
possible sites suitable for particle acceleration in order to explain the observations at very 
high energies. Here I present the results obtained from the first relativistic hydrodynamical 
simulations of jets in high-mass microquasars. Our results allow us to make estimates for 
the emission originated in diff'erent sites of the whole structure generated by the jets. These 
works represent a first step in trying to obtain a deeper understanding of the physics and 
emission processes related with jets in high-mass microquasars. 
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1. Introduction 

High-mass X-ray binaries (HMXB) are the 
parent population of high-mass microquasars 
(HMMQ), the latter being HMXB with jet ac- 
tivity. Jets of X-ray binaries (microquasars) are 
produced close to the compact object (black 
hole or neutron star) via ejection of material 
accreted from the stellar companion. Jets in 
HMMQs propagate in a medium filled by the 
stellar wind of the massive primary star of the 
system. At the scales of the binary, the stellar 
wind can significantly aflTect the dynamics of 
the jet, trigger shocks suitable for particle ac- 
celeration, and under some conditions the jet 
may be disrupted by a combination of asym- 
metric recollimation shocks and hydrodynam- 
ical instabilities (see Perucho & Bosch-Ramon 
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2008, Perucho et al. 2010). On the other hand, 
powerful jets can propagate farther, and once at 
distances from the injection point much bigger 
than the separation distance (z » <iorb, where 
z is the axial coordinate and distance to the 
compact object). Depending on the age of the 
HMXB, the jet will go through the Supernova 
remnant (SNR) created by its progenitor or, 
at later stages, the shocked stellar wind, and 
then the ISM (Bosch-Ramon et al. 201 1). After 
crossing these regions, the jet can propagate as 
far and fast as the inertia of the swept ISM per- 
mits, and eventually inflates a wide cocoon that 
pushes a shell of shocked ISM (Bordas et al. 
2009). The occurrence of collisionless shocks 
in microquasar jets can lead to efiftcient particle 
acceleration (e.g., Rieger et al. 2007, Perucho 
& Bosch-Ramon 2008, Bordas et al. 2009, 
Bosch-Ramon et al. 2011) and non-thermal 
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emission of synchrotron and inverse Compton 
origin and, possibly, from proton-proton colli- 
sions (see Bosch-Ramon & Khangulyan 2009 
and references therein). 

To understand the dynamics of these 
jets, numerical simulations of their propaga- 
tion were performed (Peter & Eichler 1995; 
Velazquez & Raga 2000; Perucho & Bosch- 
Ramon 2008, Bordas et al. 2009, Perucho et 
al. 2010, Bosch-Ramon et al. 2011) together 
with several analytical treatments (e.g., Heinz 
& Sunyaev 2002; Kaiser et al. 2004; Araudo et 
al. 2009). In this contribution, I review the re- 
sults from our simulations of microquasar jets 
at different scales. 

2. A physical picture obtained from 
the simulations 

For our study, we used a finite-difference 
code named Ratpenat, which solves the equa- 
tions of relativistic hydrodynamics in con- 
servation form using high-resolution- shock- 
capturing methods (see Perucho et al. 2010). 

2.1. Binary region 

The jet is injected supersonically (see Perucho 
& Bosch-Ramon 2008, Perucho et al. 2010) in 
the ambient and generates a bow- shock. It re- 
produces the traditional morphology of super- 
sonic outflows, i.e., the jet flow reaches the ter- 
minal (reverse) shock at the head of the jet gas 
and fills a cavity of shocked jet material that 
surrounds the beam. The cocoon is surrounded 
by the shocked ambient medium. Both jet and 
ambient shocked material mix in the contact 
discontinuity that separates the cocoon and the 
shocked ambient gas. However, inside the bi- 
nary region (z < 3 x 10^^ cm), the ambient 
consists on a stellar wind from the compan- 
ion massive star. In the case of massive stars 
(types OB), this wind can be powerful, with 
mass losses on the order of 10"^ M© yr"^ and 
typical velocities on the order of 2 x 10^ cm s"^ 
(see Perucho & Bosch-Ramon 2008 and refer- 
ences therein). In this situation, the bow-shock 
is impacted by this wind on the side where 
the companion star is with respect to the com- 
pact object at the moment of injection. On the 



other side, the bow-shock expands in the rar- 
efaction generated by the wind in the same di- 
rection. This asymmetry causes pressure differ- 
ences on both sides of the jet. As long as the jet 
is powerful (Lj > 10-^^ erg/s), and the pressure 
jump in the shock is large enough, the asym- 
metry induced in the bow- shock shape is not 
significant. Figure [T] shows the last snapshot of 
the simulation for the case of a jet with power 
Lj > 10^^ erg/s. This image shows the influ- 
ence of the wind on the shape of the bow- shock 
for weak jets and the jet helical morphology 
arising after the recollimation shock. 



The jet flow can be overpressured with re- 
spect to the surrounding cocoon and undergo 
an adiabatic expansion until it gets under- 
pressured with respect to the ambient, recol- 
limating in a strong reconfinement shock. Up 
to this shock, the large overpressure of the 
jet makes it insensitive to the asymmetry in 
pressure generated by the wind. After the re- 
confinement, the flow is decelerated and the 
aforementioned asymmetry generates helical 
Kelvin-Helmholtz instabilities that affect the 
whole body of the jet. These instabilities may 
cause the disruption of jets with power Lj < 
^q37 gj,g^g within the binary region. Disruption 
by instabilities causes mass entrainment, de- 
celeration and an increase of the jet cross- 
section. 



The winds from massive stars are known 
to be clumpy and this feature could be of 
importance regarding the jet propagation, due 
to strong mass entrainment by dense clumps 
that propagate into the jet. We have performed 
simulations of an overpressured jets in direct 
contact with the wind (i.e., we implicitly as- 
sume that the bow shock is far away and the 
cocoon/backflow has been evacuated by the 
wind). In these simulations we have included 
thermal radiative cooling terms. The jet pow- 
ers are Lj = 3x 10^^ erg/s and Lj = 10^^ erg/s. 
Our results show strong entrainment and inter- 
action with individual knots (see Fig.[2l). This, 
together with the wind thrust, causes a devi- 
ation of the jet's original direction of motion 
and efficient deceleration of the jet material. 
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Fig. 1. Cuts of logarithm of rest mass density in the plane defined by the binary star and the jet (top) and 
perpendicular to it (bottom) for a jet with power Lj - 10^^ erg/s after 1000 s. 



2.2. Outside the binary region 

Once the jet is outside the binary region, at 
z ~ 3 X 10^^ cm ^ 1 pc, the ambient that 
it will encounter will depend on the age of 
the HMXB (Bosch-Ramon et al. 2011). In the 
case of a young HMXB, with an age of 4rc ~ 
3x10^ yr, the jet would encounter the stellar 
wind shocked by the shocked SNR, which to 
its turn drives a strong shock in the ISM. In 
this case, a jet with power Lj = 3 x 10^^ erg/s 
is strongly decelerated by the shocked ISM and 
it would tipically need several thousands of 
years to cross this region and reach the ISM. 
In the case of a HMXB with an age around 
4rc -10^ yr, after crossing the shocked wind 
region, the jet is decelerated at the thin layer of 
shocked ISM. However, after a few thousand 
years, the jet head drills this layer completely 
and enters the ISM region (Fig.O. 

For an older source, with an age 4rc ~ 
10^ yr, the ISM shock velocity could become 
slower than the proper velocity of the system 
through the ISM (Vp ~ 10^ cm/s). In this case, 
a bow- shock structure would form in the di- 
rection of Vp and the shocked wind would be 
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Fig. 3. Rest mass density maps showing the inter- 
action of a jet with power Lj - 3 x 10^^ erg/s with 
the shocked stellar wind and shocked ISM. 



evacuated in the opposite direction with a sim- 
ilar velocity to its unshocked wind one v^ - 
10^ cm/s. In this case, our simulations show 
that a jet with power Lj = 3x 10^^ erg/s would 
be disrupted by the effect of the side impact of 
the shocked wind. As a result of this interac- 
tion, the jet is deviated and decelerated to the 
extent that it is not able to cross the shocked 
ISM, which acts as a wall for the jet material. 
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Logarithm of rest-mass density. Transversal to binary plane. 
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Logarithm of rest-mass density. Parallel to binary plane. 
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Fig. 2. Cuts of logarithm of rest mass density in the plane defined by the binary star and the jet (top) and 
perpendicular to it (bottom) for a jet with power Lj = 10^^ erg/s. 



2.3. At last, the ISM 

Once the jet escapes all the previously men- 
tioned interaction regions, it enters the ISM, 
where its evolution can be approximated as 
a jet propagating in a homogeneous ambient 
medium (Bordas et al. 2009). In such standard 
scenario, an overpressured jet generates a bow- 
shock in the ambient, and undergoes ones ore 
several reconfinement shocks due to its pres- 
sure difference with the surrounding cocoon. 

Bordas et al. (2009) modelled this evolu- 
tion and studied the emission from this kind 
of systems with ages between 10"^ and 10^ yrs. 
The calculations of the emission was based on 



the previous description of the system, which 
relies on both analytical and numerical work. 



3. Particle acceleration sites and 
emission 

In the previous sections, we have described 
some important processes in the evolution of 
jets in HMXBs. In all of them, strong interac- 
tions take place, which are suitable locations 
for particle acceleration: a) Bow-shock, reverse 
and reconfinement shock in the binary region 
and in the ISM, and b) strong interaction of the 
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jet with the shocked SNR and/or ISM at differ- 
ent ages of the HMXB. 

Our results put limits to the possible detec- 
tion of emission from jets in HMXBs in the 
different studied spatial and temporal scales. 
The estimates we make for the different sit- 
uations show that the radiation produced at 
recollimation shocks inside the binary region, 
~ 10^^ cm, could explain the observed high- 
energy emission from several high-mass X-ray 
binaries (HMXB) (see, e.g., Bosch-Ramon & 
Khangulyan 2009 and references therein). 

Once the jet is far outside the binary re- 
gion, and depending on its age, it could be em- 
bedded in the hot shocked wind/SNR ejecta. 
At this stage, significant non-thermal emission 
can be generated in the reverse shock and co- 
coon region, mainly radio and X-rays, and also 
very high-energy photons, via synchrotron ra- 
diation and inverse compton scattering of pho- 
tons from the companion star. On the contrary, 
as the jet propagates through the hot ambi- 
ent, there is no forward bow-like shock until 
it reaches the shocked ISM. In that stage, the 
jet is pressure confined inside the hot shocked 
wind and/or SNR and this leads to pinching 
and reconfinement. Thus, particle acceleration 
may take place in these regions. See the discus- 
sion in Bosch-Ramon et al. (2011) about the 
implications of these results for sources like 
SS433 and Cy gnus X-1. 

In the case of jets propagating in the ISM, 
Bordas et al. (2009) showed that MQs could 
be extended non-thermal emitters from radio 
to gamma-rays, with emission being gener- 
ated in the shocked ambient medium, reverse 
shock and in reconfinement shock(s) along the 
jet. Although these results strongly depend on 
the non- thermal luminosity fraction and on the 
magnetic field assumed for the calculations, 
the model is able to account for sources from 
which radiation could be detected and those 
from which it would not, mainly depending on 
jet power, age and ambient density. In general, 
the emitted flux is larger for larger jet power 
and ambient density, whereas the changes due 
to the source age are frequency and radiative- 
process dependent. 
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